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SUMMARY 


Si' 


A  theoretical  study  is  made  of  the  flight  path  of  an  uncontrolled 
but  aerodynamically  stable  satellite  as  it  spirals  down  through  the  earth's 
atmosphere  from  an  initially  circular  orbit  under  the  action  of  air  drag. 
The  earth  and  its  atmosphere  are  taken  as  spherically  symmetrical.  Par 
altitudes  above  about  125  n.  miles  a  simple  solution  is  found:  the 
velocity  of  the  satellite  is  independent  of  its  size,  shape,  weight  and 
initial  altitude  and  equal  to  the  orbital  velocity  appropriate  to  its 
current  altitude;  while  its  angle  of  descent,  in  radians  is  twice  its 
drag/weight  ratio.  Estimates  of  lifetime  are  also  made. 


1  Introduction 


Now  that  earth-satellites  are  emerging  from  the  domain  of  fancy 
into  the  colder  climate  of  reality,  the  citizens  who  live  beneath  their 
orbits  will  be  asking  sharper  questions  about  the  satellites'  probable 
behaviour  when,  under  the  prolonged  action  of  air  dra&  they  come  down 
to  earth.  A  20  lb  satellite  seems  trivial  enough,  but  the  average 
householder  might  be  somewhat  daunted  at  the  prospect  of  a  20  ton 
satellite,  red-hot,  landing  in  his  garden.  The  likelihood  of  such  a 
landing  depends  on  the  speed  and  angle  of  descent  of  the  satellite  at 
the  altitudes  where  severe  aerodynamic  heating  occurs;  and  this  in 
turn  depends  on  the  past  history  of  the  satellite,  on  how  it  behaves 
as  it  spirals  down  through  the  upper  atmosphere.  In  this  Memo,  it  is 
shown  that  at  altitudes  above  about  125  n.  miles  simple  and  general 
formulae  for  the  velocity  and  angle  of  descent  can  be  derived,  if  certain 
simplifying  assumptions ,  stated  below,  axe  accepted. 
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2  Assumptions 


The  two  main  assumptions  are: 

(i)  The  earth  is  spherical  and  non-rotating,  with  a  spherically 
symmetrical  atmosphere  and  gravitational  field;  and  the  satellite  is 
initially  in  a  circular  orbit. 


(ii)  The  satellite  is  uncontrolled  and  perfectly  stable  aerodynamically, 
jb.  \  remaining  at  zero  incidence  throughout  its  descent  path.  Thus  its  axis 
AC  cV^  ;P°ints  in  the  direction  of  motion  and  the  only  force  acting,  other  than 
'  gravity,  is  axial  drag. 
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These  assumptions  are  helpful  simplifications  in  a  first  attack 
on  the  problem  and  it  is  difficult  to  guess  how  serious  sure  the  errors 
incurred.  In  practice  the  earth's  oblateness,  or  rather  the  atmosphere's 
oblateness,  would  impose  a  periodic  variation  on  the  axial  drag.  Also, 
unless  the  satellite  was  a  sphere  its  incidence  would  not  be  zero  but 
would  vary  periodically  in  a  weathercock  oscillation  of  low  frequency; 
and  superposed  on  this  oscillation  would  be  another  resulting  from  the 
rotation  of  the  atmosphere.  Intuitively  one  would  expect  that  these 
various  disturbances  would  nearly  Impose  small  perturbations  on  the 
simplified  solution  derived  here,  though  there  is  always  the  risk  that 
the  perturbations  may  build  up  through  resonance. 

3  Analysis 

Pig. 1  shows  the  notation.  CO  is  a  fixed  line  through  the  initial 
position  of  the  satellite  and  S  is  its  position  at  any  later  time  t. 

The  two  forces  acting  on  the  satellite  are  drag  D  parallel  to  its  axis, 

2  2 

and  weight  mgR  /r  ,  where  m  is  its  mass,  g  the  acceleration  due  to 
gravity  at  the  earth's  surface ,  R  the  earth's  radius  and  r  the  distance 
CS  between  satellite  and  earth's  centre.  Let  <f>  be  the  angle  between  CS 
and  CO,  let  v  be  the  satellite's  velocity,  and  let  0  be  its  angle  of 
descent,  i.e.  the  inclination  of  its  flight  path  to  the  local  horizontal. 

Then,  resolving  tangential  and  normal  to  the  flight  path,  we  have, 
as  the  dynamical  equations  of  motion, 

»§£  =  2^  sin  8-D  (1) 

r 


{<f>  +  0)  =  cos  0  ,  (2) 


since  (0+6)  is  the  angle  between  the  direction  of  motion  and  the 
fixed  line  OB.  Equations  (l)  and  (2)  are  supplemented  by  two  kinematic 
equations  expressing  the  radial  and  transverse  velocities,  namely 


dr 

dt 


-  v  sin  0 


(3) 


V  cos  0 


(4) 


If  we  assume  that  0  is  small,  so  that  sin  6  can  be  taken  as  0 
and  cos  0  as  1 ,  and  use  equation  (4)  to  eliminate  t,  equations  (l)  - 
(3)  become 


2 

dv  _  gR  6  rD 
d0  ~  r  m 


(5) 


2  d0 
v 

C L<p 


_  v2 


(6) 
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§  ■  -rS  (7) 

Dividing  equations  (5)  and  (6)  by  equation  (7)  to  eliminate  <f> : 


dv 

dr 


2ft  d0 
v  0  ~  a 
dr 


r2  m© 

+  £ 
r2  r 


(8) 

(9) 


As  long  as  0  is  small,  the  velocity  of  the  satellite  is  likely  to 

be  near  the  orbital  speed  appropriate  to  its  altitude,  gR2/r:  for  if 

v  were  very  much  less  than  orbital  speed  the  satellite  would  promptly 
'fall  out  of  its  orbit’,  i.e.  0  would  no  longer  be  small.  So  it  is 
logical  to  make  the  substitution 


v2  =  &£  -  X 
r 

and  to  assume  that 

X  « 

r 

Eliminating  v  from  equation  (9)  by  means  of  equation  (10) 


(10) 


(11) 


£ 

=  “  (62)  to  the  first  carder  (12) 


since  X  «  gR  /r.  Differentiating  equation  (10),  we  have 


2v 


dv 

dr 


gR  dX 
r2  "  dr 


(13) 


Eliminating  v  dv/dr  from  (8)  and  (13) : 
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"by  equation  (12).  To  solve  equation  (14)  we  assume  first  that  since  0 
is  small  the  right-hand  side  of  (14)  is  almost  negligible,  so  that 


2 

2Rr 

mgR2 


D 


-.2  ,  2 
mgR  /r 


2 

weight 


(15) 


Adopting  the  usual  aerodynamic  notation  we  now  define  the  drag 
coefficient  by 


ipv  S 


06) 


where  p  is  the  density  of  the  air  around  the  satellite,  and  S  is  a 
reference  area,  taken  as  the  maximum  cross-sectional  area  of  the  satellite 
perpendicular  to  the  direction  of  motion. 


With  this  notation  equation  (15)  becomes 


if  we  denote  SC^/m  by  A.  Since  is  constant  at  high  Mach  numbers 

in  free -molecule  flow,  A  may  be  taken  as  constant.  Also  v2r/gR2  is 
nearly  constant  by  equations  (10)  and  (ll),  while  pr  varies  enormously 
with  altitude.  So,  using  (l 7),  we  may  write 


d  ,  x2 

—  (PO 

dr 


(18) 


where  f 


and  dashes  denote  differentiation  with  respect  to  r. 


Equation  (14)  may  be  re-written 


and,  on  using  (16)  and  (18)  this  becomes 


2 

or,  since  (v  r/gR  ) 


1, 


0 


pr 


£ zpt'j 


(19) 
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Also  in  equation  (12)  we  may  take 


so  that,  from  (10)  and  (12) 


=  at 

r 


rf ’A 


'A2') 


J 


(20) 
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to  the  first  order}  on  taking  v  r/gR  =  1  in  the  second  (small)  term. 

Equations  (19)  and  (20)  are  the  first  order  solutions  for  6  and  v. 

In  order  to  evaluate  f’  and  f”  we  must  specify  a  "standard 
atmosphere",  i  e.  fix  the  variation  of  p  with  r.  At  present  however 
the  air  density  at  high  altitudes  is  not  known  at  all  accurately,  and  any 
chosen  "standard  atmosphere"  will  no  doubt  prove  erroneous  when  the  true 
values  are  known.  The  values  chosen  here  for  p  between  100  and  200  n. 
miles  altitude  are  given  in  the  table  below.  With  these  values  of  p, 
the  quantity  pr  can  be  approximated  by  the  function 


pr 


^2674.  160  -  236.28Q4r  x  10"6  +  5.195r2  x  10~12 


(21) 


the  maximum  error,  for  altitudes  between  100  and  200  n.  miles,  being  about 
3/5,  as  the  Table  shows.  Values  of  rf'  and  r2f ' ' ,  found  by  differentiating 
2 

f  =  (pr)  as  given  by  (21 ),  are  also  included  in  the  Table. 


TABLE  -  The  variation  with  altitude  of  p,  pr,  etc. 


Altitude 

n.  miles 

100 

120 

140 

160 

180 

200 

p  lb/cu  ft 

471  xl  0"1 2 

-12 

103x10 

-12 

27.9x10 

- ZTS" 

8.65x10 

^12" 

3.18x10 

ITS 

1.30x10 

pr  lb/sq  ft 

1 0. 1  xl  0"3 

2. 23x1 0"3 

607x1 0-6 

1 89x1 0"6 

„  -6 
69.9x10 

28. 8x1 0-6 

Approx,  for 
pr  (21) 

1 0. 1  xl  0-3 

2. 29x1 0"3 

605x10 

188x1 0-6 

67.7x1 0~6 

28. 7x1 0"6 

-rf  ’ 

55. 6x1 0-3 

2.61x1 0"3 

1 64x1 0"6 

13.9x1 0"6 

1.56x10-6 

0. 235x1 0"6 

r2f  • ' 

31.6 

1.35 

76. 7x1  o’"3 

5. 83x1  O'3 

578x1 0"6 

75.8xlo“6 

2 

The  Table  provides  values  of  rf'  and  r  f»»  for  substitution  in 
equations  (19)  and  (20):  it  only  remains  to  choose  values  for  A.  At 
altitudes  above  about  100  n.  miles  the  mean  free  path  of  the  air  molecules 
is  greater  than  the  dimensions  of  the  satellite  and  ' free -molecule '  aero¬ 
dynamics  must  be  used.  In  these  conditions  the  drag  coefficient  of 

the  satellite,  whatever  its  shape,  may  conveniently  be  taken  as  2  as  long 
as  the  satellite's  speed  is  much  greater  than  the  mean  molecular  speed. 

So  A  =  SGp/m  =  2 S/m,  where  S  is  the  satellite  cross-sectional  area 

in  sq  feet  and  m  is  its  mass  in  lb.  Next  we  must  set  limits  to  S/m. 
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If  we  exclude  objects  like  metal-foil  'balloons',  and  also  exclude  anything 
smaller  than  the  U.  S.  20”  sphere  of  mass  20  lb,  we  may  take  0.1  as  an  upper 
limit  for  S/m.  As  an  example  at  the  other  end  of  the  scale  we  may  take  a 
10,000  lb  satellite  with  a  cross-sectional  area  of  50  sq  ft  («  8  ft  diameter), 
giving  S/m  -  0.005.  Thus  we  assume  here  that  A  lies  between  0.01 
and  0.2. 


If  the  simplest  solutions  for  v  and  0,  namely 


(22) 


and 


0  =  prA  , 


(23) 


2  2  2 

are  to  be  reasonably  accurate,  the  terms  r  f''A  /2  and  rf 'A  /2  in 

equations  (19)  and  (20)  must  be  small.  The  Table  shows  that  |rf '  |  is 

always  less  than  V300  of  r2f",  and  so  the  error  in  equation  (22)  will 
always  be  much  less  than  the  error  in  equation  (23).  If  equation  (23)  is 

2  2 

to  hold  with  say  1%  accuracy,  r  f"A  must  be  less  than  0.02,  and  using 
2 

the  values  of  r  f"  in  the  Table,  we  find  that  the  maximum  permissible 
values  for  A  are  0,025  at  100  n.  miles  altitude,  0.12  at  120  n.  miles 
altitude  and  0.51  at  14-0  n.  miles.  Thus  the  values  of  0  given  by 
equation  (23)  are  in  error  by  less  than  1%  down  to  about  127  n.  miles 
altitude  if  A  =  0.2  or  about  90  n.  miles  if  A  =  0.01.  Under  the  same 
conditions  the  values  of  v  given  by  (22)  are  in  error  by  less  than 
0. 005$.  Equation  (23)  can  be  re-written  in  the  same  form  as  (15), 


for 


since 


e  =  2  -&%g- 

weight  * 


prA 


prSCn 


?pv2scD 

mgR2/r2  * 


2 

v 


gR2/r  . 


4  Results 

At  altitudes  where  equations  (22)  and  (23)  hold,  i.e.  above  about 
90  n,  miles  altitude  if  A  is  small  (=:  0.01)  or  above  about  125  n. 
miles  if  A  is  large  (=*  0. 2),  the  speed  of  the  satellite  at  any  altitude 

is  equal  to  the  orbital  speed  appropriate  to  that  altitude,  gR2/r , 

and  the  angle  of  descent  0  (the  inclination  of  the  flight  path  to  the 
local  horizontal)  is  very  nearly  equal  to  prA  ,  or  twice  the  drag/weight 
ratio,  the  maximum  error  being  1%. 

Pig. 2  shows  how  the  speed  varies  with  altitude,  and  Pig. 3  gives  the 
variation  of  0  with  altitude  between  200  and  100  n.  miles,  accepting 
the  values  of  p  given  in  the  Table.  At  130  n.  miles  altitude  the  speed 
has  reached  25,460  ft/sec,  while  the  angle  of  descent  lies  between 
0.00066°  for  A  =  0.01  and  0.0132  for  A  =  0.2. 

The  contrast  between  Pig. 2  and  Pig. 3  is  interesting.  The  satellite's 
speed  is  independent  of  its  shape,  size  and  weight,  i.e.  is  independent  of 
A,  whereas  its  angle  of  descent  is  directly  proportional  to  A  and  hence 
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inversely  proportional  to  the  mass/area  ratio  m/S.  The  increase  in  air 
density  as  the  satellite  descends  causes  a  rapid  increase  in  the  angle  of 
descent,  which  is  proportional  to  the  density.  But  the  increase  in  drag  as 
the  satellite  descends  does  not  reduce  the  speed,  as  might  be  expected. 

On  the  contrary,  the  speed  goes  on  rising  steadily,  because  the  angle  of 
descent  always  adjusts  itself  so  that  half  the  potential  energy  accruing 
from  loss  in  altitude  is  absorbed  by  the  drag  while  the  other  half  goes  to 
increase  the  kinetic  energy.  This  increase  in  speed  as  the  satellite  descends, 
which  of  course  iirrplies  a  shorter  period  of  revolution,  does  not  go  on 
indefinitely.  The  velocity  reaches  a  maximum  (at  about  80  n.  miles  altitude 
when  A  =  0.02)  and  then  falls  rapidly. 

Once  the  angle  of  descent  0  is  known,  the  number  of  revolutions  N 
made  by  the  satellite  in  dropping  from  any  given  initial  altitude  can  easily 
be  found  by  integration  of  equation  (7)  in  the  form 


which  shows  incidentally  that  N  is  proportional  to  Vo,  i.e.  to  Va,  for 
given  r.  The  number  of  revolutions  performed  by  satellites  having  various 
values  of  A  and  an  initial  altitude  of  200  n.  miles  are  shown  in  Fig. 4, 
a  scale  of  days  also  being  included.  From  this  graph  the  number  of  revolutions 
completed  between  any  two  altitudes  (in  the  100-200  n.  miles  bracket)  can  be 
found:  e.g.  for  A  =  0.04,  N  «  750  revolutions  at  150  n.  miles,  and  N  “  810 
at  120  n.  miles,  i.e.  about  60  revolutions  are  made  in  descending  fi*om  150 
to  120  n.  miles  altitude.  It  is  clear  from  Fig. 4  that  very  few  more 
revolutions  will  be  performed  below  100  n.  miles;  so  the  time  taken  to 
drop  to  100  n.  miles  can  conveniently  be  taken  as  the  lifetime  of  the 
satellite.  Lifetime,  so  defined,  is  plotted  against  initial  altitude  for 
various  A  in  Fig. 5.  The  lifetimes  plotted  in  Fig. 5  are  in  fair  agreement 
with  those  given  in  reference  1 ,  which  were  derived  in  a  different  way. 

The  values  of  v  and  0  given  in  Figs.  2  and  3  are  those  obtained 
from  the  simplest  solutions,  equations  (22)  and  (23).  v  and  0  can  be 
found  with  similar  accuracy  down  to  an  altitude  some  7  n.  miles  lower  by 
using  equations  (20)  and  (19)  instead  of  (22)  and  (23).  The  lower  limit  of 
altitude  can  be  pushed  a  mile  or  two  further  still  if  the  second-order 
approximations  for  v  and  0  are  used.  These  are  found  to  be 


v 


f « »  +  2f 1  f  "  '  +  ff '  ‘  • 


)] 


if  certain  terms  which  can  be  shown  to  be  small  are  ignored. 
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Symbol  Unit 

Gp  =  drag  coefficient  =  D/£pv2S ,  taken  as  2. 


D 

e 

f 

g 

m 

N 

r 

R 

S 

t 

v 

A 

0 

X 

% 

P 

<f> 


^g  pdls 

2. 71 828. . . . 

(pr)2  lb2/ft^ 

acceleration  due  to  gravity  at  earth's  surface,  taken 

as  32.2  ft/sec 

mass  of  satellite  it 

number  of  revolutions 

distance  of  satellite  from  earth's  centre  ft 


earth's  radius  (taken  as  10,800/rc  n.  miles,  or 
20.9015  x  106  ft) 

cross-sectional  area  of  satellite,  in  plane 


perpendicular  to  its  axis  sq  ft 

time  sec 

velocity  of  satellite  ft/sec 

scD/m  ft2/lb 

angle  of  descent  of  satellite 

=  inclination  of  flight  path  to  local  horizontal 

gR2/**  ~  v2  ft2/sec2 

3.14159.... 

air  density  lb/ft^ 

angular  travel  of  satellite  relative  to  earth's 


centre  (see  Pig. l) 


# 
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